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Abstract 



In this paper we use the latest corrections to the Newton-Einstein 
secular rates of perihelia of some planets of the Solar System, phc- 
nomenologically estimated with the EPM2004 ephemerides by the Rus- 
sian astronomer E.V. Pitjcva. to put severe constraints on the range 
parameter A characterizing the Yukawa-like modifications of the New- 
tonian inverse-square law of gravitation. It turns out that the range 
cannot exceed about one tenth of an Astronomical Unit. We assumed 
neither equivalence principle violating effects nor spatial variations of 
a and A. This finding may have important consequences on all the 
modified theories of gravity involving Yukawa-type terms with range 
parameters much larger than the Solar System size. However, caution 
is advised since we, currently have at our disposal only the extra-rates 
of periehlia estimated by Pitjeva: if and when other groups will es- 
timate their own corrections to the secular motion of perihelia, more 
robust and firm tests may be conducted. 
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1 Introduction 

Historically, the first attempts to find deviations from the Newtonian inverse- 
square law of gravitation were performed to explain the anomalous sec- 
ular precession of Mercury's perihelion discovered by Le Verrier (1859): 
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Hall (1894) noted that he could account for Mercury's precession if the 
law of gravity, instead of falling off as 1/r 2 , actually falls of as l/r fc with 
k = 2.00000016. However, such an idea was not found to be very appeal- 
ing, since it conflicts with basic conservation laws, e.g., Gauss's Law, unless 
one also postulates a correspondingly modified metric for space. Other his- 
torical attempts to modify Newton's law of gravitation to account for the 
Mercury's perihelion behavior yielded velocity-dependent additional terms: 
for a review of them see (Gine, 2007) and references therein. Such attempts 
practically ceased after the successful explanation of the perihelion rate of 
Mercury by Einstein (1915) in terms of his general theory of relativity: an 
exception is represented by Manev (1930) who, with a 1/r 2 correction to 
the Newtonian potential, was able to reproduce the anomalous apsidal pre- 
cession of Mercury. 

It was recently realized that deviations from the Newton's inverse-square 
law could provide windows into new physics (Fuji, 1991; Fischbach et al., 
1992). Indeed, in the modern framework of the challenge of unifying gravity 
with the other three fundamental forces of Nature possible new phenomena 
could show up as deviations from the inverse-square law of gravitation. In 
general, they would occur at submillimeter length scales, but sometimes also 
at astronomical or even cosmological distances. For a review of the many 
theoretical speculations about deviations from the 1/r 2 law see (Adelberger 
et al., 2003). 

Among various parameterizations like, e.g., power-law (Fischbach et al., 
2001), a very popular, phenomenological way to account for a possible vio- 
lation of the Newtonian inverse-square law takes the form of a Yukawa-like, 
exponentially modified Newtonian potential 



r 

l + a exp ( - - 



(1) 



where G is the Newtonian gravitational constant, M is the mass of the 
central body which acts as source of the gravitational field, a and A are the 
strength 1 and the range, respectively, of the putative new interaction. The 
Yukawa correction to the Newtonian potential 

GMa ( r\ 
U Y = — exp(^--J (2) 



1 Here we will not consider composition-dependent a which would induce violations of 
the equivalence principle; for a derivation of the potential of eq. (1) in a relativistic gravity 
model obeying the equivalence principle see (Zhytnikov and Nester, 1994). 
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yields an entirely radial extra-acceleration 

GMa /■_, r\ ( r\ 
A Y = -— (l + -)exp(--). (3) 

For a review of various theoretical frameworks (braneworld models, scalar- 
tensor or scalar-tensor-vector theories of gravity, studies of topological de- 
fects) yielding a Yukawa-like, fifth force see, e.g., (Krause and Fischbach, 
2001; Bertolami and Paramos, 2005; Bertolami et al, 2006; Moffat, 2006a) 
and references therein. Among them, there are various models of modified 
gravity which predict effects at astronomical scales or even larger. For exam- 
ple, the recent Scalar- Tensor- Vector Gravity (STVG) by Moffat (2006a) in 
the intentions of his proponent would be able to comprehensively and consis- 
tently account for the observed data in the Solar System, the Galaxy, clus- 
ters of galaxies and cosmological scenarios (Moffat, 2006b). Other studies 
on long-range, Yukawa-like modifications of gravity conducted with differ- 
ent techniques on astronomical/astrophysical scales can be found in (White 
and Kochanek, 2001; Amendola and Quercellini, 2004; Sealfon et al., 2005; 
Reynaud and Jaekel, 2005; Shirata et al., 2005; Sereno and Peacock, 2006). 

The problem of finding experimental or observational constraints on the 
parameters a and A, which is crucial to exclude unviable models and achieve 
some progress in the study on those that appear feasible, is usually tackled by 
looking at what happens at a by keeping A fixed, and subsequently repeating 
the process by sampling different spatial ranges for A, without asking if this 
or that particular range for A is, in fact, really allowed: see, e.g. (Mikkelsen 
and Newman, 1977; De Rujula, 1986; Burgess and Cloutier, 1988; Fischbach 
and Talmadge, 1999; Krause and Fischbach, 2001; Bertolami and Nunes, 
2003; Sereno and Jetzer, 2006; Iorio, 2007a). 

In this paper we will test, in a purely phenomenological way, a very 
definite and widely used assumption in many modified theories of gravity, 
i.e. the hypothesis that A may assume values of the same order of magni- 
tude, or larger than the typical sizes of the planetary orbits in the Solar 
System (Moffat, 2006b). We will show that Solar System tests are, in fact, 
able to tell us something important about ranges A > 10 11 m. To this 
aim, we will, first, derive an explicit expression of the secular, i.e. averaged 
over one orbital revolution, perihelion precession induced by a Yukawa-type 
anomalous acceleration on the orbits of the Solar System planets. Then, 
we will compare our formula, obtained in the small eccentricity approxima- 
tion, with the latest estimated corrections of the perihelion 2 rates (Pitjeva, 

2 The perihelia, as the other Keplerian orbital elements, are not observable quantities: 
ranges, range-rates, right ascensions, declinations are, in fact, measured. 
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2005a) in order to see if results consistent with the tested hypothesis are 
obtained. About the methodology adopted, it is important to note that 
the corrections to the perihelion rates determined in (Pitjeva, 2005a) are 
phenomenologically estimated quantities of a global, least-square solution 
in which only Newtonian and Einsteinian dynamics was modeled: no ex- 
otic dynamical terms were included in the fit. Thus, in our opinion, such 
phenomenological corrections can genuinely be used to get information on a 
hypothetic, unmodeled force. It may be interesting to note that a somewhat 
independent test of the reliability of such a strategy can be found in (Iorio, 
2007b) in which the mass of the Kuiper Belt Objects was assessed with the 
extra-rates of perihelia by Pitjeva (2005a) by obtaining results compatible 
with other estimates from different, non-dynamical techniques. If and when 
other groups will estimate their own corrections to perihelion rates we will 
use such determinations as well in order to enforce and extend our test. 
A complementary approach which could be followed consists in repeating 
the global fit of the Solar System data by modifying the dynamical force 
models of the data reduction softwares with the addition of the investigated 
non-standard acceleration term and, accordingly, including in the set of the 
parameters to be estimated in the least-square sense the ones connected 
with the Yukawa potential as well, so to look at their mutual correlations as 
well 3 : however, such a strategy would be model-dependent and might yield 
just the outcomes desired by the experimentalist. 

2 The effects of a Yukawa-like fifth force on the 
perihelia 

To be more definite, let us suppose that a given theory, for various theoretical 
and/or observational reasons, makes use of a A quite larger than the typical 
spatial scales of the Solar System, e.g. because of a fit of a data set of a 
physical system different from it. In this case, an independent test of such 
an assumption is to check if a A with such characteristics yield, in fact, 
results compatible with the determined Solar System dynamics, within the 
associated errors. Clearly, should un-physical and/or inconsistent results be 
obtained, the considered model(s) and the related hypothesis would be ruled 
out. 

Let us, now, work out the orbital effects induced by eq. (3), treated 
as a small perturbation of the Newtonian monopole term, on the planetary 

3 According to some people, this would be the only trustable approach to the problem. 
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motions of the Solar System planets. In view of a direct comparison with 
the latest estimated extra-rates of the longitude of perihelion w, we will 
consider the secular precession of such an element. For a radial perturbing 
acceleration A r , the Gauss equation for the variation of w can be written as 



dw \f\ 



e 



2 



■Ar COS /, (4) 



dt nae 

where a is the planet's semimajor axis, e is the eccentricity, n = \J GM/a 3 
is the Keplerian mean motion and / is the true anomaly. In order to obtain 
the secular rate of w, eq. (3) must, first, be evaluated upon the unperturbed 
Keplerian ellipse, given by 

r = a(l — ecosE) (5) 

in terms of the eccentric anomaly E; then, it must be inserted into the 
right-hand-side of eq. (4) and, finally, the integral over a complete orbital 
revolution must be performed. The following formulas will be used 

r «»/ = Ss&> 

(6) 

dt = {1 - ecosE) dE, 

In the calculation, which we are going to perform by quite reasonably as- 
suming that a and A are constant and uniform over the typical spatial and 
temporal scales of Solar System bodies, the expression 

exp ^— cos Ej (7) 

appears; it prevents us from obtaining a closed form of the averaged per- 
ihelion rate because the modified Bessel functions of first kind /o,i(&e/A) 
would appear (Burgess and Cloutier, 1988). Let us assume A > ae; with 
this choice, we can safely use 

(ae \ ae 
exp y— cos E j rs 1 + — cos E. (8) 

In the small eccentricity approximation we, thus, get 

aVGMa / a\ 
^^^exp(--), (9) 

up to terms of order 0{e 2 ). Expressions analogous to eq. (9) can be found 
in (Burgess and Cloutier, 1988; Talmadge et al., 1988; Reynaud and Jaekel, 
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2005), in which quantities proportional to the perihelion advance after one 
orbital revolution were worked out 4 , and in (Sereno and Jetzer, 2006) where 
the the perihelion secular rate was calculated up to 0(e ). 



3 Constraining the range and the strength of a 
Yukawa-like fifth force with planetary perihelia 

The formula of eq. (9) is very useful because it allows us to get important 
information on the size of A. Indeed, let us write down eq. (9) for a pair of 
planets, say A and B, and take their ratio: by assuming that both a and A 
do not vary with distance we get 



zu 



ZU 



(B) 




cxp 



«B - a A 
A 



(10) 



Note that the ratio of the rates of perihelia due to a Yuakawa-like interaction 
is independent of a not only when the approximated expression of eq. (9) is 
used, but also when the general expression with the Bessel function (Burgess 
and Cloutier, 1988) is adopted. By defining 



n 



zu 



(A) 



(B) 



and 

e(A) = 

it is possible to construct 




OA / Ob - a A 

exp - 

a B V A 



T(A) = n-9(A); 



(11) 

(12) 
(13) 



if, for a given range of values of A > ae|A/B |T I turns out to be incompatible 
with zero, within the errors, that range for A must be discarded. Note that 
our analysis is independent of a, assumed to be nonzero, of course. The 
uncertainty in T can be conservatively assessed as 



<5T(A) < 5U + 56(A), 



with 



5u < mi 



+ 



zu 



(A) I 



37(B)! 



ZU 



(14) 



(15) 



4 They are 2-kzu/u (Burgess and Cloutier, 1988; Talmadge et al., 1988) and vj/n (Rey- 
naud and Jaekel, 2005). 
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Table 1: Estimated semimajor axes a, in AU (1 AU= 1.49597870691 x 10 11 
m) (Pitjeva, 2005b), and phenomenologically estimated corrections to the 
Newtonian-Einsteinian perihelion rates, in arcseconds per century (" cy _1 ), 
of Mercury, the Earth and Mars (Pitjeva, 2005a). Also the associated errors 
are quoted: they are in m for a (Pitjeva, 2005b) and in " cy" 1 for w (Pitjeva, 
2005a). For the semimajor axes they are the formal, statistical ones, while 
for the perihelia they are realistic in the sense that they were obtained from 
comparison of many different solutions with different sets of parameters and 
observations (Pitjeva, private communication 2005). The results presented 
in the text do not change if 5a are re-scaled by a factor 10 in order to get 
more realistic uncertainties. 



Planet 


a (AU) 


5a (m) 


w (" cy" 1 ) 


5m (" cy'") 


Mercury 


0.38709893 


0.105 


-0.0036 


0.0050 


Earth 


1.00000011 


0.146 


-0.0002 


0.0004 


Mars 


1.52366231 


0.657 


0.0001 


0.0005 



5Q(X) < 6(A) 



2aA 



5aA + 



1 1 

2ae A 



5a B 



(16) 



The linear sum of the individual errors in eq. (15) accounts for the exist- 
ing correlations among the estimated perihelia corrections, which reach a 
maximum of about 20% for Mercury and the Earth (Pitjeva, private com- 
munication 2005). For A=Earth and B=Mercury, Table 1 and Figure 1 tell 
us that A ~ aMer is not allowed at about 1.5 — a level; for larger heliocentric 
distances the constraints are quite tighter, exceeding the 3 — a level. Al- 
though the inspected range for A ends at Mars in Figure 1 , it turns out that 
larger values, far beyond the Solar System boundaries, are ruled out as well 
at about 8 — a level. Our results are unaffected by re-scaling by a factor 10 
the formal errors in the semimajor axes. 

After having discovered that A cannot exceed the semimajor axis of Mer- 
cury, let us now further constrain it. From eq. (10) it can be obtained 



A 



a B - a A 



In 



1A 



(17) 



It turns out that the major sources of error are the estimated extra-rates of 
perihelia through their ratio II, so that 



5X < 



a B - oa 



nin^ 



1A 



<5n. 



(18) 
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5= 3.5 



1.5 




0.5 1 1.5 

MAU) 



Figure 1: |T|/JT from the data of the Earth and Mercury over a range 
dMer < A < aMar- Values of A > aMar are ruled out at an even larger 
number of a. 
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For A=Earth, B=Mercury we have 



A = 0.182 ± 0.183 AU, (19) 

which is marginally compatible with zero. Note that the use of eq. (9), from 
which eq. (10) and eq. (17) come, can be a posteriori justified because for 
the Earth and Mercury the obtained value for A yields oe/A < 1. 

The result of eq. (19) allows us to constrain a as well. Indeed, in the 
case of Mars we have 

y = 0.78, (20) 

so that the approximation of eq. (8), and the formula of eq. (9) based on 
it, hold. Thus, from eq. (19) and the values of Table 1 for Mars we get 

a= . exp - = 2 x 10" 10 . (21) 







a 




2 - 








A 



VGMm, va 
The uncertainty can be evaluated as 

5\ + mSw ] = 1.3 x 10~ 9 . (22) 
\w\ J 

Also in this case, a is compatible with zero. Such constrains on a are less 
tight than those obtained, e.g., in (Iorio, 2007a; Sereno and Jetzer, 2006), 
but the authors of such works made use of values of A which the present 
analysis has ruled out. 



4 Conclusions 

In this paper we put on the test the hypothesis that modifications of the 
Newtonian inverse-square law, parameterized in terms of a Yukawa-like cor- 
rection, can occur over astronomical scales by using the corrections to the 
Newtonian-Einsteinian secular rates of the perihelia of Mercury and the 
Earth phenomenologically estimated, in the least-square sense, with the 
EPM2004 ephemerides by Pitjeva (2005a). 

By taking their ratio we found that the range parameter A of a Yukawa- 
like fifth force cannot exceed about 0.18 AU. The determined extra-precession 
of the perihelion of Mars yielded an upper bound on a of 10~ 9 , which is com- 
patible with other estimates obtained with other approaches. The values 
obtained for both A and a are compatible with zero; moreover, the results 
presented here are left unaffected by re-scaling the uncertainties in the esti- 
mated Keplerian orbital elements by a factor 10 in order to evaluate them 
more realistically. 
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If and when corrections to the secular rates of perihelia will be estimated 
by other teams of astronomers, more complete and extensive tests could be 
performed. 

Another approach which could be followed consists in introducing an ad- 
hoc Yukawa-type term in the dynamical force models of the ephemerides- 
generating routines and repeating the global fit of the whole Solar System 
data set by estimating, among other things, also the parameters in terms of 
which the new force is expressed. 
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